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Abstract 

We present analytic estimates for the energy levels of N electrons 
(N = 2 — 5) in a two-dimensional parabolic quantum dot. A magnetic 
field is applied perpendicularly to the confinement plane. The rele- 
vant scaled energy is shown to be a smooth function of the parameter 
P= (effective Rydberg/effective dot energy) 1 / 6 . Two-point Pade ap- 
proximants are obtained from the series expansions of the energy near 
the oscillator (J3 — > 0) and Wigner (/? — > oo) limits. The approximants 
are expected to work with an error not greater than 2.5% in the entire 
interval < (3 < oo. 
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1 Introduction 



The problem of N electrons in quantum dots and magnetic fields has been 
widely considered in recent years [1]. 

Concerning the energy levels of N electrons in model two-dimensional 
parabolic quantum dots, the actual magnitude to compute is a scaled energy 
which depends only on one parameter, (3— (effective Rydberg/effective dot 
energy) 1 / 6 . When (3 — * (very high magnetic fields, for example), we may 
use perturbation theory to compute the energy levels. In the opposite limit, 
(3 — > oo, a strong coupling expansion may be used to obtain the energy. 
The idea of the present paper is to construct two-point Pade approximants 
interpolating from (3 — to (3 — > oo. We show results for 2, 3, 4, and 5 
electrons. 

To our knowledge, there are only a few alternative analytical ways to 
obtain the energy of certain levels in the entire interval < (3 < oo. Semi- 
classical [2] and 1/| J|-expansions [3, 4], both working for states with high 
angular momentum, J, are available. Besides these expansions, there is also 
an idea of improving the perturbative series by using the asymptotics at 
(3 — > oo [5]. We understand the present paper as a useful complement to the 
results of [2 - 5] . We will see, for example, that it is very simple to find Pade 
approximants for states with \ J\— 0, 1, 2, for which the methods of [2, 4] do 
not work. 

2 Two-point Pade approximants 

The construction of approximants follows the idea of paper [6], in which the 
hydrogenic energy levels in a magnetic field were obtained. 

Let us consider the expansions of the scaled energy (to be computed in 
the next sections) when (3 — > and (3 — > oo, 



s 



e 



J2hf3 k + 0(f3 s+1 ), 



(1) 



fc=0 



e 



/3->oo 




(2) 
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Many of the coefficients entering (1, 2) are zero. For example, all the bk with 
k 7^ mod 3, a±, a 3 , etc. 

A two-point Pade approximant is a rational function, 

reproducing the expansions (f) and (2). go may be fixed to one. The asymp- 
totics when (5 — ■> oo forces L to be equal to K + 2. Equating number of 
coefficients in the expansions (1) and (2) to number of unknowns in (3), we 
obtain, 

s + t = 2K+l. (4) 

At a given K, there is a set of possible pairs s, t. We will attach indexes 
s and t to the approximant, P s t . Among the possible P Sit , the best one takes 
nearly the same number of terms in the expansions (1) and (2) [6], i.e. s ~ t. 
For example, at K = 3 the best one, which coefficients are computed almost 
trivially, is P^. 

Let us consider the equations fulfilled by the pu and the qu- Equating (3) 
to (1) and (2), we obtain, 



Pk = Qoh + Qih-i + • • • + QkK, < A; < Min(s, K + 2), (5) 

Pk = QKdK+2-k + <iK~l^K+l-k + ■ ■ ■ + Qk-2d0, 

K + 2 — Mm(K + 2, t) < k < K + 2. (6) 

Eliminating the overlapping pk, i.e. k in the interval Max(0, s + 1 — K) < 
k < Min(s, K + 2), a system of linear equations for the coefficients qt arise. 
We assume that s is in the interval K — 1 < s < K + 2, such that this system 
contains K equations. Let us write explicitly, for example, the approximant 
P 4 3 . The coefficients 61, b 2 , 64, a\ and a 3 are assumed to be zero. The result 
is the following, 



hP 3 a2 L 1 + qiP 



where q 2 = a /(b - a 2 ), qi = a q 2 /h, q 3 = (a qi - b 3 )/(b - a 2 ). (7) 
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Formula (7), or similar expressions for higher approximants (see Appendix 
1), is to be used throughout the paper. Let us indicate the way to compute 
the bk and the a k . 



3 The coefficients bk and 

Let us consider the two dimensional motion of N electrons in a parabolic 
quantum dot of energy Houq- A magnetic field is applied normaly to the plane 
of motion. The hamiltonian governing the internal motion (center of mass 
motion is excluded) is written in dimensionless variables as 

H , j 9^c q fa , 

m = h + 2n J+ 2n s - (8) 



where u c is the cyclotronic frequency, Q = ^juq+u^/A is the effective dot 
frequency, J is the total (internal) angular momentum (along the z axis), S z 
is the z-projection of the total spin, g is the effective giromagnetic factor, 
and 



~ h W + Pk d P J h U + plJ del 
+ h pU 9o k de k+1 + n - i ^ [pi pi +1 ) do k 



Notice that h depends only on one parameter, f3 3 = ^ pe A / ' k 2 Ti 2 / 'VhQ, where 
p is the effective electron mass, and k - the relative dielectric constant. The 
coordinates entering h are the moduli of the Jacobi vectors, 



Pk = J— . k = l,...,N-l, (10) 




and the angles between p k and Pk+i, denoted by 9 k . The dimensionless 
reduced masses are p k — k/(k + 1). 
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The eigenfunctions of H are e iJ ~ij), where S accounts for global rotations, 
and ip are the eigenfunctions of h. The eigenvalues of H are trivially ob- 
tained from the eigenvalues of h, which will be called e. We will obtain Pade 
approximants to e. 

In the (3 — > (oscillator) limit, perturbation theory may be applied to 
obtain e. The resulting series is the following 

e\^ = b + b 3 p 3 + b 6 [3 6 + ..., (11) 

where b = N — 1 + |J| + 2n, n is the total number of oscillator quanta, 6 3 = 

01 J2k<i r fc« 1 |0); e tc- 111 systems with more than two electrons, sometimes 
egenerate perturbation theory shall be used to compute b 3 , b e , etc. 

On the other hand, when /3 — > oo a strong coupling expansion may be 

applied. Distances are scaled according to p — > Pp. The hamiltonian h takes 

the form 



h 2 1 

? = 4 fc ?/' + S^ + (^ 



+ ^4 E , 2 




In the leading approximation, /3 — > oo, we shall minimize the classical 
potential energy entering the r.h.s. of (12). It is found that the electrons sit 
at the corners of a regular polygon. This is a few-body version of the Wigner 
solid. The classical energy becomes a function of one variable, pi, 

jj — M | BinTT/iV 1 

8sin 2 7r/iV + Pl ^|sin0 fc ,/2|' 1 J 

where is the angle between particles and /, measured from the cm. The 
minimization of U leads to the equilibrium value of pi, p w . The equilibrium 
values of the other coordinates are obtained from geometric considerations. 
The energy in this approximation is given by U(p\o). Corrections to this 
value are obtained by writting p k = pko + Zk/P, k = 1, . . ., N — 1, 9 k = 
&ko + ZN-i + k/(3, k = 1, . . . , iV — 2, and expanding 
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h/f3 2 = U( Pw ) + h 2 /f3 2 + h 3 /f3 3 + ..., (14) 
i> = ^o + ^i/(3 + ^ 2 /(3 2 + ..., (15) 
e/(3 2 = a + a 2 /f3 2 + a 3 /f3 3 + .... (16) 

The Schrodinger equation is split into the set of uncoupled equations 



a = U(p w ), (17) 
h 2 ip = «2^o, etc. (18) 

The hamiltonian h 2 describes harmonic oscillations around the Wigner struc- 
ture. The expression for a 2 is thus 

2N-3 

a 2 = E WfcK + 1/2), (19) 
k=i 

where the Uk are the normal frequencies. Higher corrections are obtained by 
considering h 3 , /14, etc as perturbations to h 2 . 

Below, we present results for 2, 3, 4, and 5 electrons. 



4 Two electrons 

As mentioned above, the eigenfunctions of H are written as e iJ= -?/>(Pi)- Under 
a permutation of particles, S changes by ir, and does not change. Thus, 
even \J\ are related to antisymmetric spin functions, S = 0, and odd \J\ are 
related to states with spin 5 = 1. 

Let (f>k be the eigenfunctions of h at f3 — 0, 

fc = C fcP i J| 4 J| (p?/2)e^/ 4 , (20) 

where C fe = (2' J I(A;+ | JDI/A;!)" 1 / 2 , and L^" 7 ' are generalized Laguerre polyno- 
mials. The corresponding eigenvalues are 



6o = 2ib + |J| + l. (21) 

We take a fixed k as the unperturbed level, let us say, n\. Higher coefficients 
of the expansion are computed from 
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h = <n 1 \ — \n 1 >, (22) 



1 „ | < m|l/pi|A; > | 2 

66 = <--*> ■ etc - (23) 

We show in Tab. 1 the coefficients 6 3 and 6 6 for a set of two-electron 
states. 

Let us consider now the opposite limit, (3 — > oo. The equilibrium value 
of pi (scaled) is p w = 2 1 / 3 . The coefficient a is thus a = 3/2 4 / 3 . Then, we 
write pi = pio + z 1 /(3, and expand the hamiltonian. The results are 



to = -7^2 +~A, (24) 



ft? 2 ' 4' 

h (-l)V (-l) fc ' 3 , fc -3 9 , J 2 (-l) k -\k - 3) 

Mfc - fc+l ^1 fc-2 Z l ~nT~ "I ~fc^2 Z l ' K - ^ Z0 J 

PIO r 10 ^1 r 10 



Notice that the hk with odd contain an odd number of creation and 
anihilation operators. Thus, all the a& with odd k will be zero. Computation 
of the matrix elements of hk is a trivial task. Finally, we obtain 



a 2 = wi(ni + l/2), (26) 
a 4 = 2~ 2/3 {(n 2 + n 1 + 7/6)/6 + J 2 - 1/4 1 , etc (27) 

where uj\ = y/3, i.e. the classical result [7]. Note that we have used the same 
number, ni, to label the state at f3 — and (3 — > oo. As level crossings can 
not occur, the first state at /3 = should be the first when (3 — > oo, and so 
on. Note also that J appears for the first time in a 4 . 

The coefficients a 2 and a 4 are also shown in Tab. 1 for a set of states. 
From this Table and formulae contained in Appendix 1, we may construct 
Pade approximants for the energy levels. 

We show in Fig. 1 the three curves e|^_„ = b + b 3 j3 3 + b 6 j3 6 , e]^^ = 
a (3 2 + a 2 + a 4 //3 2 , and P 6i5 (/3) for the first state with quantum numbers 
| J | = 3, n\ — 0. This is the typical behaviour of the approximants. 
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In Fig. 2, we compare the approximants in the sequence Pr+i,k- The 
same state as in Fig. 1 is studied. The relative differences between consecu- 
tive approximants are shown. The (3 axis is compressed to (0,1). We see that 
the maximal difference reduces by a factor of two when K is increased by 
one. These results suggest the Pg,5 approximant (the highest we computed) 
to be accurate to about three parts in 10 3 or better. Notice that the maxima 
are reached at f3 ~ 2, i.e. in the region where the approximants jump from 
the weak-coupling to the strong-coupling regimes (see Fig. 1). 

Similar results are obtained for the other states in Tab. 1. When J = 0, 
however, relative differences between approximants increase up to 2%, and 
some approximants can not be used as they exhibit a pole. So, the best of 
our Pade estimates are expected to work with an error not greater than 2% 
in the worst situation. 

We show in Fig. 3 the convergence of the sequence Pr+i,k at particular 
(3 values, at which exact solutions are available [8]. It can be easily shown, 
for example, that 

1> = pf (l+Pi/yVl + l) e-rf (28) 

are eigenf unctions of h at (3 — (2| J| + l) 1 ^ 6 , with eigenvalues e = | J\ + 2. 
The corresponding n\ are zero. When < \J\ < 3, (3 is in the interval 
1 < (3 < 1.38, well outside the exactly solvable limits. The relative error of 
the P 6i5 approximant is lower than 2 % at J = 0, and less than 0.2 % at 
\J\ = 3. 



5 Three electrons 

First, let us consider the computation of the coefficients up to a^. The 
Wigner configuration is a triangle with side pio = 3 1 / 3 (scaled). The leading 
approximation to the energy is 

3 5/3 

a = — . (29) 

Then, we write p 1 = p i0 + (u + v)/(V2(3), P2 = Pio + (u - v)/(V2(3), 
Q\ = tc/2 + y^2z/(pioP). Expanding the hamiltonian h, we arrive to 
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= - ho + Tb + ^i )+-y + W + z*), (30) 



/i 3 = = — i -32-- + 32w-^- - 8u a - Quv z - hv 6 - 6uz 2 




(31) 

. d n . T d 

h 4 = —< — +u- o M ^ + 7 M + v a 2lJv TT 

ou 2 oz 2 4 ov oz 

_3 2 92 

2 V dz 2 

-^uvz 2 - t|ttuV + ^-z 4 \ , etc. (32) 
16 128 256 J v ; 

u, v, and z are normal coordinates, u corresponds to a symmetric oscil- 
lation (breathing) with frequency u\ = \/3, whereas v and z correspond to 
a mixed oscillation of the Wigner structure, with frequency cj 2 = ^3/2. In 
the harmonic approximation, the spatial wave function is written as 

i> = Xn 1 { U )Xn 2 {v)Xn :) {zy J *, (33) 

where the x are oscillator functions, i.e. Hermite polynomials multiplied 
by gaussians. A function ip that can be antisymmetrized is related to a 
spin-polarized state (S = 3/2) of three electrons. On the other hand, a 
mixed-symmetry is related to a spin one half state. 

Let us consider the lowest state with a given J, i.e. that one with numbers 
nx — n-i — nj, — 0. This state may be spatially antisymmetrized only when 
J = 3k, with k - an integer. The argument goes along the lines sketched 
in [9]: a cyclical permutation of electrons in the triangle, leaving the wave 
function invariant, is equivalent to a 27r/3 rotation, which multiplies it by 
e tJ2n / 3 , thus e* J27r / 3 = 1. The excitations with n 2 = n 3 = correspond also to 
antisymmetric states. The n 2 + n 3 = 1 states correspond to mixed-symmetry 
doublets, etc. On the other hand, when J ^ 3k the lowest state and the 
excitations with n 2 = n 3 = have mixed symmetry. An antisymmetric and 
a mixed-symmetry state appear at n 2 + n 3 = 1, etc. We will restrict the 
analysis to the lowest state and the first excitations. 

The coefficient a 2 is thus given by 
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o 2 = y/3(m + 1/2) + yj3/2(m + n 2 + 1). (34) 

Higer corrections are computed from perturbation theory around h 2 . By 
the same reason as for two particles, a 3 and a 5 are equal to zero. The next 
nonzero coefficient is 



a 4 = < n 1 ,n2,n3\h 4: \n 1 ,n2,n 3 > 

+ < ni, n 2 , nslh^ki, fa, fa >< fa, fa, fa\h 3 \n 1 ,n 2 , n 3 > 

(fei,fe2,fc3)^(ni,n2,n 3 ) V^Tll - fa) + ^/2(rii + U 2 - fa - fa) 

(35) 

The computation of matrix elements entering (35) is trivial, leading to 

a 4 = A ) (-2 + 9V2 + 72 J 2 + (12 + 18^)^1 + 12n 2 

+ (36 + 9^ + 18V2ni)(n 2 + n 3 ) + 228n 2 n 3 - 78(n 2 2 + n 2 3 )} . (36) 

Let us stress that we expanded the hamiltonian around a structure with 
Q\ = 7r/2. There is an equivalent configuration with 6\ = —n/2. In the 
expansion (16), we have neglected tunneling between the two equivalent con- 
figurations. The same comment holds for systems with more than three elec- 
trons. Notice also that the second local minimum of U, the linear structure 
(second "Lagrange" solution), is at a distance 0.36/9 2 above the lowest state. 
We can disregard any effect comming from this structure when (3 2 3> a 2 /0.36. 

The computation of the coefficient 63 requires the wave functions at (3 = 0. 
They may be explicitly written with the help of the indications of paper [10]. 
We will compute 63 for a set of states with | J\—0, 1, 2, and 3. In all of these 
states, the corresponding quantum numbers at (3 — > 00 can be specified. 

For example, when J = 0, the first antisymmetric (A) state at f3 — 0, 
which starts from 6 — 4 (it will be labelled (4,A)), goes to the first A at 
(3 — > 00, i.e. (ni,n 2 ,n 3 ) = (0,0,0). The first mixed state, (4,M), goes to 
the first doublet with n\ = 0, n 2 + n 3 = l.The second antisymmetric state, 
(6,A), goes to the (1,0,0). At j J| = 1, the (3,M) goes to the (0,0,0), the (5,A) 
goes to a state with n\ = 0, n 2 + n 3 = 1. At | J\ = 2, the (4,M) goes to the 
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(0,0,0), and the (6, A) - to a state with ri\ = 0, n 2 + = 1. At \ J\ = 3, the 
(5, A) goes to the (0,0,0), and the (5,M) to a state with ri\ = 0, n 2 + 113 = 1. 

Let us write explicitly the needed wave functions at (5 — (up to nor- 
malisations). In the case of mixed symmetry, only one representative of the 
doublet is given. 

J = 

4>a,a = PiPasinflx e-^)/ 4 , (37) 
04,m = (Pi - Pi) e-^'\ (38) 
<f>e,A = P1P2 sin 0! (pi + pj-8) e -^ + ®l\ (39) 

J = 1 

03,m = Pi e^/ 2 e -^l\ (40) 
05,A = {(-2 Pl p 2 + p\) e-^" - Pl pl e 3 ^/ 2 } e-^/\ (41) 

J = 2 

04,m = Pip 2 e-("? + ^)/ 4 , (42) 
6iA = p 1 p 2 sin^ 1 ( p 2 e-^+p 2 e^) e^^ 4 . (43) 

J = 3 



5 ,a = {p^/ 2 -3 Pl p 2 e^ 2 } e-^Dl\ (44) 
5 ,m = {p 3 e- 3 ^/ 2 + p lP 2 e ^/ 2 } e-rt + M*. (45) 

The computation of 63 is thus a trivial task. We groupped all the results 
together in Table 2. With these coefficients, we construct the approximants 
^3,2, -P4,3, and P 5j4 . 

The convergence analysis of the Pade sequence Pk+i,k is shown in Fig. 
4 for the lowest antisymmetric state with |J| =3. As in the two-electron 
problem, convergence is strong, suggesting the P^^ interpolant to be accurate 
to about 6 parts in 10 3 or better in the entire interval < (3 < 00. Similar 
results are obtained for the lowest states with |J| =2, and 1. The excited 
states and the states with J = show a slower convergence. The relative 
error is estimated as 2%. 
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6 Four electrons 



First, let us compute the coefficients a and a 2 . As in previous cases, a± = 
a 3 = 0. 

The equilibrium configuration at j3 — > oo is a square with side pio = 
(2 + \/2/2) 1//3 (scaled). The equilibrium values of the other coordinates are 
P20 = \ZV3 pio, P30 = Spio/v^, 6>io = arctan2 , 9 20 = 37r/4 - 9 W . a Q is given 
by 

a = 3(2 + V2/2) 1/3 . (46) 

Expanding around the equilibrium geometry, we obtain the quadratic 
hamiltonian, 



112 \dz\ + dz\ + dzl) \pj + p|oJ dz\ 

( 1 1 \ 2 9 9 T/ 

" + + + (47) 

where 



V 2 = 0.543374^ + 0.157283^2 + 0.520305^ + 0.182039^3 

+ 0.377402^3 + 0.460201z| + 0.171862^4 + 0.192148^4 

+ 0.517308^ + 0.110937^5 + 0.085931z 2 z 5 + 0.562904z 4 z 5 

+ 0.8536134 (48) 

The normal frequencies are easily found, resulting in 



ux = 1.04969, 
uj 2 = 1.13911, 
lu 3 = 1.36227, 
uj a = 1.46201, 

u 5 = 1.75640, (49) 
and, thus, the coefficient a 2 is given by 
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5 

a 2 = 5>*(n fc + l/2). (50) 

k=l 

Notice that the frequency corresponding to the breathing mode, uj 5 , is 
close to the classical value >/3, but it does not exactly coincide with \f?>. 

As in the A = 3 problem, the lowest state with a given \J\, i.e. n& = 0, 
k = 1, ... ,5, can be spatially antisymmetrised only when \J\ takes certain 
values. The allowed values are \J\ = 2,6,10,.... These are polarised spin 
states, i.e. with total spin S = 2. The excitations of the lu 5 mode have the 
symmetry of the ground state. 

Let us stress that, as N increases, the number of equivalent configurations 
(geometries with the same ao) increases, and the energetic distance to other 
local minima of the classical energy decreases. Thus, tunelling effects become 
more and more important. 

In what follows, we restrict the analysis to the lowest spin-polarised state, 
that is | J| = 2, and the n& = 0. It seems that it is the lowest state of h in the 
sector with S = 2 at any f3. Indeed, as (3 — > 0, this state goes to an oscillator 
state with 6 = 7, i.e. to the (7,A) in the terminology used above. There is a 
second (7,A) with J = 0. However, according to Hund's rule, this state has 
a higher energy at f3 <C 1. When (3 — > oo, the J = state is also higher in 
energy because J = is not compatible with = 0, thus excitation quanta 
are needed. 

Let us compute the coefficient 63 for the \J\ = 2 state. The wave function 
at f3 = is given in Appendix 2. Calculations may be carried out analytically, 
leading to 

< 4>7 a\ — \4>7 A > 91 , 

= 6 V7 f p?: = A^- (si) 

<<M07,a> 64 1 ; 

Once computed ao, 02, bo, and 63, we may construct the approximants 
and P^3, which are the main results of this section. The relative error 
between P 3i2 and is shown in Fig. 5, suggesting that P^ may estimate 
the energy with an error not greater than 2% at intermediate values of (3. 

Pade approximants to other levels may be constructed in the same way. 
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7 Five electrons 



We follow the same programme as in the N = 4 problem, i.e. computation 
of do, a.2-, bo, and 63 and, from them, construction of the approximants P^2 
and P^s. 

The equilibrium configuration at (3 — > 00 is a pentagon with side pio = 
1.30766 (scaled). The equilibrium values of the other coordinates are P20/P10 — 
1.44177, pso/pio = 1.53244, p 40 /Pio = 1.34500, 1O = 0.865925 rad, # 20 = 
1.74428 rad, # 30 = -1.03941 rad. a is given by 

a = 9.28013 . (52) 

Expanding around the equilibrium configuration, we obtain the hamilto- 
nian h 2l 



ho 



where 



2 




P20 P30 ^6 0^7 



\/ 2 = 0.624037^- 0.007988^2 + 0.514889^ + 0.107324^3 

+ 0.231789^3 + 0.438330^ + 0.073587^4 + 0.245901^4 

+ 0.381949^4 + 0.373258*1 + 0.315672^^5 - 0.296995z 3 z 5 

+ 0.451706^ + 0.313498^6 - 0.217440^6 - 0.286484z 3 z 6 

+ 0.326409^-26 + 0.765563^^6 + 1.43839^ - 0.126796^7 

+ 0.119294^3^7 + 0.169244^^7 + 1.37849^^7 + 1.46079^. (54) 

The normal frequencies following from the eigenvalue problem for /i 2 are 



wi = 0.727516, 
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Uj = 



0J2 = 



u 3 = 



0U4 = 



u 5 = 



Uq = 



0.804763, 

1.37185, 

1.52505, 

1.68899, 

1.73205, 

1.74375, 



(55) 



and the coefficient 02 is given by 
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«2 = 



uj k {n k + 1/2). 



(56) 



k=l 



The allowed values of |J| for antisymmetric states with n k = are \J\ = 

0,5,10, In these states the total spin is S = 5/2. In what follows, we 

consider the lowest state in this sequence, i.e. J = 0. This state goes to a 
(10,A) state as (3 — > 0. It is not, however, the lowest of all antisymmetric 
levels at (5 1 because a second (10, A) state with \J\ = 2 minimises the 
Coulomb repulsion. 

The wave function of the J = state is given in Appendix 2. Calculations 
may also be performed analytically, leading to 



We show in Fig. 6 the relative difference between P 3i2 and P 4i3 . This 
difference is not greater than 2.5%. 

Pade approximants to other levels may be constructed in the same way. 

8 Concluding Remarks 

We have studied systems of 2 - 5 electrons in a two-dimensional parabolic 
quantum dot. The potentials involved in this problem (harmonic plus Coulomb 
repulsion) are very gentle, and lead to a smooth dependence of the energy e 
on the coupling constant (3. This fact is graphically represented in Fig. 1, 
where it is seen that the "regions of convergence" of the perturbative and 
the strong-coupling series "intersect" with each other. 




149 



(57) 
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The degrees of homogeneity of the potentials are also important factors 
towards the smoothness of e. They lead to expansion series containing only 
powers of (3 3 at (3 — > 0, and inverse powers of (3 2 at (3 — > oo. From the calcu- 
lational point of view, it means that with the help of trivial computations, 
not beyond first order perturbation theory, we may construct approximants 
up to P4 3, that is, a quotient between a 5th order and a 3rd order polyno- 
mial in (3. These approximants are exact in both f3 — > and (3 — > 00 limits, 
leading to errors not greater than 2.5 % in the small transition region where 
they jump from one expansion to the other. The accuracy may be improved 
by computing higher approximants, as shown for two and three electrons. 
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10 Appendix 1. The approximants used in 
the paper 

Pz,M = b + a [3 2 {l - (1 + h(3/a + a (3 2 /(b - a 2 ))- 1 } . (58) 
P 4 3 is given in Eq. (7). 

p (n s , , (60 - a2)q^ 

bAP) " ° l+ gi P + ...+q 4 p 4 l + qi p+... + q ^ 

+ «o/? 2 ( 1 57} , (59) 



where, 



1 

92 = a q 3 /b 3 , q x = — {b 3 + (b - a 2 )q 3 }, 

a 

a J . . b - a 2 \ 
Q4 = — < -1 + — 7 93 > , (60) 



a 4 [ &3 
and q 3 is determined from the equation 



b 3 qi + a q 2 - (bo - a 2 )q A = 0. (61) 



p (R , , hp (&o-a 2 ){g4/3 4 + g 5 /35} 

6,5 ^ j b0+ l + qi (3 + ... + q 5 ^ l + qi (3 + ... + qb (3s 



where, 



<?3 = ^{^4 - &e}> <?2 = 7—^ — {a + a 4 g 4 }, 
63 6 - a 2 

Qi = 7-{aog2 - (bo ~ a 2 )q 4 }, 

<?5 = — {h + (po - a 2 )q 3 - a gi}, (63) 
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and g 4 is found from 

- hq 2 + ao<?3 - (bo - a 2 )q 5 = 0. (64) 

11 Appendix 2. Wave functions for 4 and 5 
electrons at (3 = 

We show the explicit form of the functions 7j a, N = 4, and 0io,a, N — 5. 
The f} are measured from the c. m., and (a x b) z denotes the z component 
(normal to the plane) of the vectorial product. 

<j>7,A = {[p3 ■ n/p 3 +t(p3 X fi) 2 /p 3 ] 2 [f 2 X f 3 + f 3 X f 4 + f 4 X F 2 ] 2 

- [p 3 • r 2 /p 3 + «(p 3 x r 2 ) z /p 3 ] 2 [r 3 X f 4 + f 4 X fi + fi x f 3 ] z 
+ [P3 • r 3 /p 3 + «(p 3 x r 3 ) 2 /p 3 ] 2 [f 4 x fi + fi x r 2 + f 2 x f 4 ] 2 

r 2 [fx xf 2 +f 2 xf 3 + f 3 X f^} e -(P?+P|+Pl)/4 e 2i(ei+2e 2 )/3 > (65) 
e -(p?+^+^+^)/4 . f2 ) (fi x f^ [ fg x f4 + f4 x f g + fs x f^ 

(f 2 • r 3 )(f 2 x f 3 ) z [f 4 x f 5 + f 5 x fi + fi x f A ] z 

(f 3 • r 4 )(f 3 x r 4 ) 2 [f 5 x rl + fi x f 2 + f 2 x f 5 ] 2 

(r 4 • f 5 )(f 4 x f 5 ) 2 [fi x f 2 + f 2 x f 3 + f 3 x fi] 2 

(f 5 • n)(f 5 x fi) 2 [f 2 x r 3 + f 3 x r 4 + f 4 x f 2 ] 2 

(fi ■ f 3 )(fi X f 3 ) 2 [f 4 xf5 + f5Xf 2 + f 2 X f 4 ] 2 

(f 2 • r 4 )(f 2 x f 4 ) 2 [f 5 x fi + fi x f 3 + f 3 x f 5 ] 2 
(f 3 • f 5 )(f 3 x f 5 ) 2 [fi x f 2 + f 2 x f 4 + f 4 x fi] 2 

(fi • f 4 )(fi x f 4 ) 2 [f 5 x f 2 + f 2 x f 3 + f 3 x f 5 ] 2 

(f 2 ■ f 5 )(f 2 x f 5 ) 2 [fi xr 3 + r 3 xr A + r A x fi] J . (66) 



+ 
+ 
+ 
+ 



+ 
+ 
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Figure and Table Captions 



Tab. 1. The first nonzero coefficients bk and for a set of two-electron 

states 

Tab. 2. The first nonzero coefficients bk and for certain three-electron 

states 

Fig. 1. Weak and strong-coupling expansions (dashed lines), and the P 6i5 
approximant (solid line) for two electrons in a state with | J| = 3, n\ — 

Fig. 2. Relative differences between consecutive approximants in the se- 
quence Pr+i,k, with K ranging from 2 to 5. The state with \J\ = 3, 
n\ = is shown 

Fig. 3. Convergence of the Pade sequence Pr+i,k at f3 = (2| J| + l) 1 / 6 . Two 
electrons in states with ri\ — are considered. Solid lines represent the 
exact solutions found in [8] 

Fig. 4. The same as in Fig. 2, now for the first antisymmetric state of three 
electrons with | J\ = 3. K ranges from 2 to 4 

Fig. 5. Relative difference between P 4i3 and P 3i2 for the lowest antisymmet- 
ric state of four electrons (| J\ = 2) 

Fig. 6. The same as in Fig. 5, now for the lowest antisymmetric state of 
five electrons with angular momentum J = 
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Table 1 



J , n 


&o 


h 


K 


a 


a 2 


CZ4 


0, 
0, 1 
0, 2 


1 

3 
5 


1.253314 
0.939986 
0.802904 


-0.345655 
-0.087406 
-0.041161 


1.190551 
1.190551 
1.190551 


0.866025 
2.598076 
4.330127 


-0.034998 

0.174989 

0.594963 


1, 
1, 1 
1, 2 


2 
4 
6 


0.626657 
0.548325 
0.499367 


-0.032762 
-0.016339 
-0.010112 


1.190551 
1.190551 
1.190551 


0.866025 
2.598076 
4.330127 


0.594963 
0.804950 
1.224924 


2, 
2, 1 
2, 2 


3 
5 
7 


0.469993 
0.430827 
0.402676 


-0.011153 
-0.007016 
-0.004920 


1.190551 
1.190551 
1.190551 


0.866025 
2.598076 
4.330127 


2.484844 
2.694831 
3.114805 


3, 
3, 1 
3, 2 


4 
6 
8 


0.391661 
0.367182 
0.348211 


-0.005533 
-0.003911 
-0.002950 


1.190551 
1.190551 
1.190551 


0.866025 
2.598076 
4.330127 


5.634647 
5.844634 
6.264608 



Table 2 



J , Symmetry 


bo 


h 


a 


a 2 


a 4 


0, A 


4 


1.87997 


3.12013 


2.09077 


0.0358156 


0, M 


4 


2.58496 


3.12013 


3.31552 


-0.0619104 


0, A 


6 


1.76247 


3.12013 


3.82282 


0.200926 


1, M 


3 


2.81996 


3.12013 


2.09077 


0.276191 


1, A 


5 


1.82122 


3.12013 


3.31552 


0.178464 


2, M 


4 


2.23247 


3.12013 


2.09077 


0.997315 


2, A 


6 


1.58623 


3.12013 


3.31552 


0.899589 


3, A 


5 


1.70372 


3.12013 


2.09077 


2.19919 


3, M 


5 


2.26184 


3.12013 


3.31552 


2.10146 
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